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^ ; Abstract 
Oh. . . 

I , The like-sign dimuon charge asymmetry of the B meson, which was reported in the D0 Collaboration, 
Q-i' 

1^ ■ is studied in the SU{5) SUSY GUT model with S4, flavor symmetry. Additional CP violating effects from 

■ 

the squark sector are discussed in Bs — Bg mixing process. The predicted like-sign charge asymmetry is in 
' the 2a range of the combined result of D0 and CDF measurements. Since the SUSY contributions in the 
quark sector affect to the lepton sector because of the SU{fy) GUT relation, two predictions are given in 
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\ the leptonic processes: (i) both BR(/i — t- 67) and the electron EDM are close to the present upper bound, 
10 • (ii) the decay ratios of r decays, r — )■ /i7 and r — )• 67, are related to each other via the Cabibbo angle Ac: 
■ BR(r — )■ e7)/BR(T — )• /U7) ~ A^. These are testable at future experiments. 



X 



; PACS numbers: ll.30.Hv, 12.60.Jv, 13.20.Hc, 14.40.Nd 



* E-mail address: ishimori@musc.sc.niigata-u.ac.jp 

^ E-mail address: kajiyama@muse.sc.niigata-u.ac.jp 

^ E-mail address: shimizu@muse.sc.niigata-u.ac.jp 

§ E-mail address: tanimoto@muse.sc.niigata-u.ac.jp 



I. INTRODUCTION 



The CP violation in the K and Bd mesons has been well explained within the framework of the 
standard model (SM). There is one phase, which is a unique source of the CP violation, so called 
Kobayashi-Maskawa (KM) phase in the quark sector with three families. Until now, the KM 
phase has successfully described all data related with the CP violation of K and Bd systems. 

However, there could be new sources of the CP violation if the SM is extended to the super- 
symmetric (SUSY) models. The CP violating phases appear in soft scalar mass matrices. These 
contribute to flavor changing neutral currents (FCNC) with the CP violation. Therefore, we should 
examine carefully CP violating phenomena in the quark sector. 

The Tevatron experiments have searched possible effect of the CP violation in the B meson system 
[21, 0] . Recently, the D0 Collaboration reported the interesting result of the like-sign dimuon charge 
asymmetry Ali(D0) = —(9.57 ±2.51 ± 1.46) x 10"'^ [s*]. This result is larger than the SM prediction 
yl^;(SM) = (— 2.3lQg) X 10~^ [3] at the 3.2a level, which indicates an anomalous CP violating phase 
arising in the Bg meson mixing. 

Actually, new physics have been discussed to explain the anomalous CP violation in several 
approaches. As a possibility, new physics contribute to decay width of the Bs meson j^- 12|. Another 
possibility is to assume new physics does not give additional contribution to the decay width but 

This typical model is the general SUSY model with gluino-mediated 
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22| . Relevant mass insertion (MI) parameters and/or squark 



the Bs — Bs mixing 13 1- 
flavor and CP violation 

mass spectrum can explain the anomalous CP violation in the Bg system. Since the squark flavor 
mixing is restricted in K and Bd meson systems, the systematic analyses are necessary to clarify the 
possible effect of squarks. 

In this paper, we study the flavor and CP violation within the framework of the non-Abelian 
discrete symmetry 27|] of quark and lepton flavors with SUSY. Then, the flavor symmetry controls 
the squark and slepton mass matrices as well as the quark and lepton ones. For example, the predicted 
squark mass matrices reflect structures of the quark mass matrices. Therefore, squark mass matrices 
provide us an important test for the flavor symmetry. 

The non-Abelian discrete symmetry of flavors has been studied intensively in the quark and lepton 
sectors. Actually, the recent neutrino data analyses 28|- 31| indicate the tri-bimaximal mixing {32! 
- 35] , which has been at flrst understood based on the non-Abelian flnite group jaSfiol- Until 
now, much progress has been made in the theoretical and phenomenological analysis of ^4 flavor 
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model |41l |- 

An attractive candidate of t he flavor s ymmetry is the 5*4 group, which was successful to explain 
both quark and lepton mixing 112| - 149| . Especially, 5*4 flavor mo dels to un ify quarks and leptons 



have been proposed in the framework of the SU (5) SVSY GUT |ll8Ml2l| . SO{W) SUSY GUT 
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124| . and the Pati-Salam SUSY GUT |l25l . Il26l |. These unified models seem to explain both 



mixing of quarks and leptons. 

Some of us have studied 5*4 flavor model 119| . which gives the proper quark flavor mixing angles 
as well as the tri-bimaximal mixing of neutrino flavors. Especially, the Cabibbo angle is predicted 
to be 15° due to 6*4 Clebsch-Gordan coefficients. Including the next-to-leading corrections of the 5*4 
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symmetry, the predicted Cabibbo angle is completely consistent with the observed one. 

We give the squark mass matrices in our 6*4 flavor model by considering the gravity mediation 
within the framework of the supergravity theory. We estimate the SUSY breaking in the squark 
mass matrices by taking account of the next-to-leading S4 invariant mass operators as well as the 
slepton mass matrices. Then, we can predict the CP violation in the Bg meson taking account of the 
constraints of the CP violation of K and mesons. We also discuss the squark effect on 6 — )■ 57 
decay and the chromo-electric dipole moment (cEDM) . 

Since our model is based on SU (5) S USY GUT, we can predict the lepton flavor violation (LFV), 
e.g., /i — )■ 67 and r — > /i7 processes |147] . In particular, the r — )■ /i7 decay ratio reflects the magnitude 
of the CP violation of the Bg meson. 

This paper is organized as follows: In section 2, we discuss the possibility of new physics in the 
framework of the CP violation of the neutral B system. In section 3, we present briefly our S4 
flavor model of quarks and leptons in SU{5) SUSY GUT, and present the squark and slepton mass 
matrices. In section 4, we discuss numerically the CP violation of the Bg meson with constraints of 
flavor and CP violations of K and B^ mesons. We also discuss the EDM of the electron, cEDM of 
strange quark and LFV. Section 5 is devoted to the summary. In appendices, we present relevant 
formulae in order to estimate the flavor violation and the CP violation. 



II. Bs - Bs MIXING 

In this section, we briefly discuss the theory and experimental results of the CP violation of the 
neutral B meson system. The effective Hamiltonian T-Ll^iq = d,s) of Bg — Bg system is given in 
terms of the dispersive (absorptive) part M'^{r'^) as 

ni^ = - '-T^, (1) 

where the off-diagonal elements M^, and Tf, are responsible for the B^ — Bq oscillations. The light (L) 
and heavy {H) physical eigenstates with mass M^^s^ and the decay width ^^h) obtained 

by diagonalizing the effective Hamiltonian "H^g. The mass and decay width difference between B^ 
and Bjj are related to the elements of Tilg as 

AM, ^Mjj-Ml = 2\Ml,l AP, ^ P^ - P^ = 2\V\,\ cos0„ 0, = ^ig{-MlJV\^, (2) 

where we have used APg ^ AMg. 

The "wrong-sign" charge asymmetry d^^ of Bg ^ fi~ X decay is defined as 

^ TiB, /.^X) - TjB, ^-X) ^ (VIA _ £k sinc^ (3) 



The like-sign dimuon charge asymmetry A^i is defined and related with a^i as 150 1 



JV++ — M — 

= /V++ , = (0-506 ± 0.043)af, + (0.494 ± 0.043)a:„ (4) 

b b 



3 



where is the number of events of bb /i^/i^X. 



±,,± ■ 



The SM prediction of A^i IS given as 

A^SM) = (-2.3^°!) X 10-^ 

which is calculated from ^ 

af;(SM) = (-4.81};^) x 10"^ ali{SM) = (2.06 ± 0.57) x lO'l 
Recently, the D0 collaboration reported A^i with 6.1 fb~^ data set as 

Ali{D0) = -(9.57 ± 2.51 ± 1.46) x 10'^ 



(5) 



(6) 



(7) 



which shows 3.2 a deviation from the SM prediction of Eq.([5]). On the other hand, the result by the 
CDF collaboration with 1.6 fb"^ data \Z] /l^;(CDF) = (8.0 ± 9.0 ± 6.8) x lO'^ is consistent with the 
SM prediction while it has large errors. Combining these measurements, one can obtain 



^^;(CDF + D0) = -(8.5 ± 2.8) x 10'^, 



which is still 3 cr away from the SM prediction. 



152 



as a:,(D0) 



(8) 



'1.7± 



The D0 Collaboration have performed the direct measurement of a^i 
9Ati^) X 10^'^, which is consistent with the SM predictio n because of its large errors. However, if 
one use the present experimental value of af; Il53 , 
find that 0,E3 



1541 1 . af;(exp) = -(4.7 ±4.6) x 10"^, one can 



all = -0.0146 ± 0.0075, 



(9) 



is required to obtain A^iiDQ). The central value of the required |a^;| is about three orders of 
magnitude larger than the SM prediction a^;(SM). Combining all results, one can obtain the average 
value 



a^j(average) ~ -(12.7 ±5.0) x 10" 



(10) 



which is still 2.5 a away from the SM prediction a^j(SM). Therefore, if the D0 result is confirmed, 
it is a promising hint of new physics (NP) beyond the SM. 

The contribution of NP to the dispersive part of the Hamiltonian is parameterized as 



where the SM contribution Mff^^ is given by 



q,SM 
12 



12772 



9' 



12) 



^ Recently, the SM predictions are updated jl5l| by the same authors. However in this paper, we use the widely- 
accepted results of Ref. 
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Input 




Input 






(231 ± 3 ± 15) MeV 


Bs{mb) 


0.841 ±0.013 ±0.020 


fSs/fBd 


1.209 ±0.007 ±0.023 




1.01 ±0.01 ±0.03 


VB 


0.8393 ± 0.0034 


Soixt) 


2.35 


Ms 


5.3663 ± 0.0006 GeV 




17.77 ±0.10 ±0.07 ps-i 


Md 


5.27917 ± 0.00029 GeV 


AAffP 


0.507 ±0.005 ps-i 




(5.1 ± 1.3) X 10-3 GeV 


nisimb) 


0.085 ±0.017 GeV 




4.248 ± 0.051 GeV 


A-pSAI 


(0.096 ± 0.039) ps-i 


4>d,SM 


i-W.ltl-l) X 10-2 


4>s,SM 


(±7.4l°|) X 10-3 


|-pS,SM| ,1 T^,rS,SM 1 
I-"- 12 I/K'^-'i2 I 


(4.97 ±0.94) X 10-3 




i52.6tlli) X 10-4 



TABLE I: Parameters of the neutral B meson mixing and quark masses 0, [155] 



with parameters hsted in Table 1. 

Using these parameters, the mass difference of Bg meson, AMg, is given by 



AMr, = AM, 



SM 



AM^^' \A 



(13) 



Since the SM contribution to the absorptive part is dominated by tree- level decay h ccs, one 



can set = T'i2^\ In this case, the wrong-sign charge asymmetry a^i is written as 



"si 



\-pq,SM, 
I-"- 12 1 


sin (0g,sM 


+ 0Ag) 






A. 





(14) 



Taking the experimental value AMJ^^p into account, one finds that |As| is strongly constrained in 
the region \As\ = 0.92 ± 0.32 j^. Therefore, unphysical condition sm(^s,SM + 0As) = —2.56 ± 1.16 



is re quired to obtain 1 a range of the charge asymmetry (See also 156| ). Also as discussed in Ref 
157 1, by using the SM prediction of r^2* experimental values of AM^^s, they found in model- 
independent way that the like-sign charge asymmetry is bounded as — < 3.16 x 10-3, where the 
CP violation Sj/^Ks Sjj^^ are also taken into account. 

Now we discuss how to avoid this unphysical condition to obtain large charge asymmetry. As 
the first possibility, one can consider the NP contributions to PJj, which come from additional 
contributions to deca y pro cesses h — )■ ccs, t^t~ s, etc. By using the D0 a nd C DF experimental data 



of Bs — )■ J /iljcj) decay |l58| . one can subtract AP^ and [3. 



AP. = ±(0.154^n?o) V^'\ fill"" = (0.39;°:}« 



-0,/2 as (l53| 
or (1.1 



-0.14\ 

-o.isJ 



(15) 



where the sign of AP^ is still undetermined, and positive (negative) sign corresponds to the first 
(second) region of pi^'^'^ . Comparing them with the SM predictions, one finds that there still can 



^ See also Ref. 



154j for recent results. 
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exist additional contributions to from NP. This possibility has been studied in several models^ 
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In Ref ■ jl3| ■ while there are no NP contributions to their model, they employed the ex- 

perimental value of ATg of Eq.f lT5|) since there must exist theoretica l uii certainties. In the {hg,ag) 
parametrization of NP, the best fit values of [hg, as) are obtained as 162 | 



{hs,as) ~ (0.5,120°), (1.8,100°), (16) 

by taking AMg, ATg, S^k and Sj/^^ into account, with varying \T\2\ in the r ang e — 0.3ps~^. In 
that paper, one can read that the region of < is favored as seen in Refs. 14j-|l6|. 

However in ordinary SUSY models, gluino-squark box diagrams do not give additional con- 
tributions to since such diagrams do not generate additional decay modes of bottom quark. 
Therefore as the other possibility, constraint for AMg is relaxed in Refs. [itI . EsI . In those papers, 
they consider models that NP does not give additional contributions to r^2) but to MI2. They 



A,< AMd,,/AMf/ < 14 



7 | a nd t he UTfit [l9|] allowed region 



18|. See also Refs. |20|422|, y, |25l ll57| for other possibilities. 



take a conservative constraint Ci 
0.776 < AMrf,,/AM|f < 1.162 

In this paper, we consider the NP contribution to Bg — Bg mixing by gluino-squark box diagrams 
in a SU (5) SUSY GUT model with 5*4 flavor symmetry. As shall be discussed in the next section, the 
soft SUSY breaking terms and related MI parameters {6^^)ij{A, B = L, R) obey 5*4 flavor symmetry. 



In such SUSY models, there are no new contributions to |l4Hl8l. l2ll |22| . While the SUSY 
contributions to Bg — Bg mixing are induced by (5^^)23, it is constrained by 6 —t- 57 decay. Since 
the other MI parameters of down- type squark sector are related to {6^^)23 due to 6*4 symmetry, K 
and Bd meson mixing, which are affected by {6^^)12 and respectively, should also be taken 

into account. The CP violation in Bg meson system is related to cEDM of the strange quark d'^ as 
well. Moreover, the leptonic processes such as r — > /i7 affected by {6^^)23 should also be taken into 
account due to SU{5) GUT relation. 

Taking the above processes into account, we assume the following conditions in our numerical 
calculation: (i) the meson mass differences satisfy 



AMdg 
0-6 < , „ ,exp < 1-4, 
AM,7 



K,SUSY I 
12 I 

AM^^P 



< 1, 



llmM 



K,SUSY< 
12 I 



v^AM^'P 



<eK = 2.2 X 10" 



:i7) 



(ii) cEDM of the strange quark is constrained by the neutron EDM as 163l . Il64 | 



|e(if I < 1.0 X lO"^'' ecm, 
(iii) the NP contribution to the branching ratio (BR) of 6 — > 57 is constrained as 

BR(6 ^ 57)^^ < 1.0 X 10"^ 



(19) 



However, the NP contributions to will be strongly constrained by the lifetime ratio tbJtb^- We would like to 
thank A. Lenz for pointing out this point. 
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(Ti,r2) 




-^2,-^3) 












SU{5) 


10 


10 


5 


1 


1 


5 5 45 1 






S4 




2 


1 


3 


2 


r 


1111 






Z4 




— i 


-1 




1 


1 


11-11 






U{1)fn 


1 








1 





-1 






(xi,X2) ix 


i,X4) ix 


5,X6,X7) (X8,X9,X10) (Xll,Xl2,Xl3) Xl4 


SU{5) 




1 


1 


1 


1 




1 


1 




5*4 




2 


2 


3' 


3 




3 


1 








— i 


1 


—i 


-1 




i 


i 


U{1)fn 




-1 


-2 













-1 



TABLE II: Assignments of SU{5), S4, Z4, and U{l)F]\f representations. 
Whi l e th e upper bounds of LFV decay processes £i — )■ £,7 and the electron EDM are given by 
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166 | 



BR(/i ^ 67) < 1.2 X 10"^\ BR(r ^ /i7) < 4.4 X 10"^ (20) 
BR(r ^ 67) < 3.3 x 10-^ \ede\ < 1.6 x lO'^^ 6cm, (21) 

we do not take these bounds into account in the numerical calculation below. Instead, in the allowed 
parameter region of our model which can explain the like-sign charge asymmetry, we will obtain the 
predictions for LFV processes. 

We perform numerical analysis in the section IV after introducing the S4 flavor model in the next 
section. 



III. THE S4 FLAVOR MODEL 



We briefly review S4 flavor model of quarks and leptons, which was proposed in |ll9l |. As the 
model is based on SU{5) SUSY GUT, it gives sfermion mass matrices as well as quark and lepton 
mass matrices. 



A. CKM Mixing 

In the SU{5) GUT, matter fields are unified into 10 and 5-dimensional representations as 10 C 
{Q, u"^, e^) and 5 C (rf'^, L). Three generations of 5, which are denoted by Fj {i = 1, 2, 3), are assigned 
to 3 of 5*4. On the other hand, the third generation of the 10- dimensional representation, T3, is 
assigned to 1 of S4, and the first and second generations of 10, (Ti,T2), are assigned to 2 of 5*4, 
respectively. Right-handed neutrinos, which are SU{5) gauge singlets, are also assigned to 2 for 
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the first and second generations, {N^,Np, and 1' for tlie tliird one, A^^. Tlie 5-dimensional, 5- 
dimensional, and 45-dimensional Higgs of SU{5), H^, ifg, and H^^ are assigned to 1 of 5*4. In order 
to obtain desired mass matrices, we introduce SU{5) gauge singlets Xj? so called flavons, which couple 
to quarks and leptons. 

T he Z4 . symmetry is added to obtain relevant couplings. Further, the Froggatt-Nielsen mecha- 



167l | is introduced to get the natural hierarchy among quark and lepton masses, as an additional 



nism 

U{1)fn flavor symmetry, where denotes the Froggatt-Nielsen flavon. The particle assignments of 
SU{5), S4, Z4^, and U{1)fn are presented in Table UTI 

The couplings of fiavons with fermions are restricted as follows. At the leading order, {xs, Xa) 
are coupled with the right-handed Majorana neutrino sector, (xs, xe, Xr) are coupled with the Dirac 
neutrino sector, (xsjXqjXio) and (xn, X125 X13) are coupled with the charged lepton and down-type 
quark sectors. At the next-to- leading order, (xi,X2) are coupled with the up-type quark sector, and 
Xi4 contributes to the charged lepton and down-type quark sectors, and then the mass ratio of the 
electron and down quark is reproduced properly. 

Our model predicts the quark mixing as well as the tri-bimaximal mixing of leptons. Especially, 
the Cabibbo angle is predicted to be 15° at the leading order. The model is consistent with the 
observed CKM mixing angles and CP violation as well as the non-vanishing Ues of the neutrino 
flavor mixing. 

Let us write down the superpotential respecting 5*4, Z4 and U{1)fn symmetries in terms of the 
S'4 cutoff scale A, and the U{1)fn cutoff scale A. In our calculation, both cutoff scales are taken as 
the GUT scale which is around lO^^GeV. The SU{5) invariant superpotential of the Yukawa sector 
up to the linear terms of Xi = 1, ■ ' ' ; 13) is given as 

w = y^TuT^) ® T3 ® (xi, X2) ® H,/A + y^T^ ^n®H, 
+ y^{N:,Np®{N:,Np^QyA 
+ y^iN:, AT^) ® (A^:, ® (X3, X4) + MAT; ® N!; 
+ y^{N:, N;) ® (Fi, F2, F3) (X5, X6, Xr) ^H,^ 0/(AA) 
+ y^N'; ® (Fi, F2, F3) ® (xs, Xe, Xi) » H^/^ 
+ yi{F,, F2, F3) ® (Ti, T2) ® (X8, X9, Xio) ® ^45 ® e/(AA) 
+ y2{F,, F2, F3) ® T3 ® (xii, X12, X13) ® H-,/A, (22) 

where y", 1/2 , yi , y^, , y2, yi, and j/2 are Yukawa couplings of order one, and M is the right-handed 
Majorana mass, which is taken to be lO^^GeV. 

In order to predict the desired quark and lepton mass matrices, we require vacuum alignments for 
the vacuum expe ctation values (VEV's) of fiavons. According to the potential analysis, which was 
presented in 119|, we have conditions of VEV's to realize the potential minimum {V = 0) as follows: 

iXi, X2) = (1, 1), (xs, Xi) = (0, 1), (X5, X6, Xr) = (1, 1, 1), (X8, X9, Xw) = (0, 1, 0), 
(Xii,Xi2,Xi3) = (0,0, 1), Xm = -— X?, (23) 

where these magnitudes are given in arbitrary units. Hereafter, we suppose these gauge-singlet scalars 
develop VEV's by denoting (xi) = a^A. 
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Denoting Higgs doublets as hu and hd, we take VEV's of following scalars by 

{hu) = Vu: {hd) = Vd, (/I45) = ^^45, (0) = 0, 



(24) 



which are supposed to be real. We define X = 9/Ato describe the Proggatt-Nielsen mechanism ^. 

Now, we can write down quark and lepton mass matrices by using the 5*4 multiplication rule in 
Appendix A. The down- type quark mass matrix at the leading order is given as 

/ \ 



yiXagV45/V2 yiXagv^s/VQ 

y2ai3Vd ) 



(25) 



Then, we have 



M\Md 







ll/iAag 



iTjI^iAagp 
||^iAag|2 









(26) 



where we denote y\Vd — yiV45. This matrix can be diagonalized by the orthogonal matrix uf^ as 

/ cos 60° sin 60° o\ 



U, 



(0) 



-sin 60° cos 60° 
V 



(27) 



The down-type quark masses are given as 



ml 



0, ml 



- \ |2 2 z 



I |2 2 2 

ml disVd 



(28) 



The down quark mass vanishes, however tiny masses appear at the next-to-leading order. 
The relevant superpotential of down sector at the next-to-leading order is given as 

= VaATi, T2) (g) (Fi, F2, F3) (g) (xi, X2) <E) (Xii, X12, X13) «) 
+ yA,(Ti, Ta) ® (i^i, F2, F3) (8) (X5, X6, Xt) ^ Xu 
+ yA^Ti, T2) (8 (Fi, F2, Fa) ® (xi, X2) ® (Xs, Xe, X7) ® ^^45/A' 
+ |/A,(Ti, T2) ® (Fi, F2, F3) ® (xii, X12, X13) ® Xl4 ® i^45/A' 
+ l/Ae^a ® (Fi, F2, F3) ® (X5, X6, X7) ® (X8, X9, Xio) ® ® /A' 

+ |/AyF3(8) (Fi,F2,F3) (g) (X8,X9,Xlo) ® (Xll,Xl2,Xl3) ®i?45®/A^ , 

which gives the correction terms in the down-type quark mass matrix. 
The down-type quark mass matrix including the next-to-leading order is 



(29) 



Md 



( eu 

\/3ms 1 



621 



V 



ei3 



£31 

2 I £32 
623 rrib + 633 J 



ei2 ^ + 622 



(30) 



^ Notice that this A is not related to the Cabibbo angle Ac in our model. 
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where the exphcit forms of e^j's are given in Appendix B, and and rrib are given in Eq. (l28i) . This 
mass matrix can be diagonahzed by 



(0)rr{i) 



(31) 



where mixing matrices for left-hand U^'' and for right-hand are estimated as 



U. 



(1) 



V 



ad 
^12 



nd nd nd 

od I nd nd 



nd 
-t>23 



"23 
nd nd i 
C^12D'i3 -L , 



V, 



a\ 

d 



(32) 



1 



where a denotes the typical value of the square root of the relevant sum of aja^'s as discussed in the 
next subsection. We neglect CP violating phases then mixing angles 6f2, Of 3, are given as 

0d o f"^] =0 (0.05) , 9f. = (—] = (0.005) , ei, = (—] = (0.005) . (33) 
\msj ynibj \mbj 

On the other hand, the superpotential of up sector at the next-to-leading order is 

Awu = ylSTi.T^) ® {T,,T2) ® ixi,X2) ® ixi,X2)®H,/A^ 
+ vl^iT^, T2) ® (Ti, T2) ® xi4 ® Xm ® HrJA^ 
+ yin ® T3 ® (X8, X9, Xio) ® (X8, X9, Xio) ® 



(34) 



Then the mass matrix becomes 



6*^14 



V 



(35) 



This symmetric mass matrix is diagonahzed by the unitary matrix Uu as 

/ cos45° sin 45° o\ /l o\ 



\ 



sin 45° cos 45° 
1/ 



Tt r 
\0 -rc rtj 



(36) 



where Tc = ^/rric/ [rric + m^) and = \/mt/ {rric + rrit). 
Therefore, the CKM matrix l^can be written as^ 



v = ul 



/l o\ 

e-'P 
\o ij 



^d ^d 



(37) 



^ The renormalization group effect for the CKM matrix is smaU so that the matrix given in the text can be regarded 
as the one at the electroweak scale. 
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where the phase p is an arbitrary parameter originating from complex Yukawa couphngs. 

At the leading order, the Cabibbo angle is derived as 60° — 45° = 15° and it can be naturally 
fitted to the observed value by including the next-to-leading order as follows: 

Vus ~ cos 15° + sin 15°. (38) 

The Vcb and Vub mixing elements are expressed as 

Vud ^ cos 15° - {ef^ + ef^e^^) sin 15°, 

Vcb ^ -nei^e'P sin 15° + r^^^gC^^ cos 15° - , (39) 
Vub ^ cos 15° + sin 15°, 

which are consistent with observed values. 

We can also predict mass matrices of the charged leptons and neutrinos, which give the tri- 
bimaximal mixing of leptons. Details are shown in Appendix C. 

Since mass eigenvalues of quarks and leptons are give in terms of Oj = (xi)/A, we can estimate Oj 
by putting the observed quark and lepton masses. These are given as 



as = ag = aio = an = ai2 = 0, ai = a2 



(yf A)2(m^3 - m^Jm^2^ V^m^^M 
a4 — ——2 , a5 — as — a/ 





2 




_ 

'2 ^2 





V^\yMvd y2Vd 



= /-I _ I, — , ai3 = , (40) 



where masses of quarks and leptons are given at the GUT scale, and the light neutrino masses mjy^ ^ 3 
are given in the Appendix C. Hereafter, we take A = 0.1 in our calculations. 



B. Squark and slepton mass matrices 

Here we study SUSY breaking terms in the framework of 6*4 x Z4 x U{1)fn to derive squark and 
slepton mass matrices. We consider the gravity mediation within the framework of the supergravity 
theory. We assume that non-vanishing F-terms of gauge and flavor singlet (moduh) fields Z and 
gauge singlet fields (i = 1, • • • , 14) contribute to the SUSY breaking. Their F-components are 
written as 

F^" = -e^K^"^ (^diW + , (41) 

where Mp is the Planck mass, W is the superpotential, K denotes the Kahler potential, Kjj denotes 
second derivatives by fields, i.e. Kf j = djdjK and K^'' is its inverse. Here the fields $fc correspond 
to the moduli fields Z and gauge singlet fields Xi- The VEV's of are estimated as = 

0(777.3/2), where 7773/2 denotes the gravitino mass, which is obtained as 7773/2 = {e'-^^'^'''^^W / M"^) . 
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First, let us study soft scalar m asses . Within the framework of the supergravity theory, the soft 
scalar mass squared is obtained as |l68[ 



LM 



(42) 



The invariance under the S^x Z4 x U{1)fn flavor symmetry as well as the gauge invariance requires 
the following form of the Kahler potential as 



(43) 



1=1,2,3 



i=l,2 



at the lowest level, where Z^^\^) and Z'^^^''^^^{^) are arbitrary functions of the singlet fields By use 
of Eq. fH2|) with the Kahler potential in Eq. f H3|) . we obtain the following matrix form of soft scalar 
masses squared for 5 5'^ and 10 lO'^ combinations, which are denoted as and m^, respectively: 



/ ml 



F)ij 



\ 



ml, 
\ mj, J 



( 



[m. 



TJij 



m. 



T{1) 





V 



m. 





2 

T{1) 





m. 



\ 



2 , 
T(2) / 



(44) 



That is, three right-handed down-type squark and left-handed slepton masses are degenerate, and 
first two generations of other sectors are degenerate. These predictions would be obvious because 
the three generations of the F C {df^, L) fields form a triplet of 5*4, and the T C {Q, u'^, e'^) fields form 
a doublet and a singlet of S'4. These predictions hold exactly before S^x Z^x U{1)fn is broken, but 
its breaking gives next-to-leading terms in the scalar mass matrices. 

Next, we study effects due to S4 x Z4 x U{1)fn breaking by Xi- That is, we estimate corrections 
to the Kahler potential including Xi- Since (Ti, T2) are assigned to 2 and its conjugate representation 
is itself 2. Similarly, (^1,^2,^3) are assigned to 3 and its conjugation is 3. Therefore, for the -^12,3 
fields, higher dimensional terms are given as 

AKp = zZ]mF,, F2, Fs) ® {F^,F^,FS) ® {x^,X^+l) ® (x?,Xm)/A' 



i=l,3 



+ Y F2, F3) ® (F^, F^, F^) ® (x., X^+l, X^+2) ® ix-, X-+2)/A' 



j=5,8,ll 



+ zZ\^m, F2, Fa) ® (F^, F^, F^) ® xu ® Xu/^' 
+ zZ\^m, F2, F3) ® (F^, F^, F^) ® ® eVAl 



(45) 



For example, higher dimensional terms including (xi,X2) and (xs^XeiX?) are explicitly written as 



AK 



A2 

A2 



(F2F3* + F3F2* + F1F3* + F3F; + FiF; + F2F*) 



(46) 
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When we take into account corrections from all XiX^ to the Kahler potential, the soft scalar masses 
squared for the -Fi 2 3 fields have the following corrections, 



F)ij 



where kp is a parameter of order one, and a%j.{k = 1, 2, 3) are linear combinations of aja/s. 
For the Ti 2,3 fields, higher dimensional terms are given as 



1=1.3 



+ ^Z] ^2) ® (Tr, ® (Xm X.+i, X^+2) ® (X- , X-+2)/A' 

i=5,8,ll 

+ z£)($)(Ti, T2) ® (Tl^ T^^) ® xi4 ® Xm/A' 



+ zfj($)(Ti, T2) ® ® (xi, X2)/A' + Z'l'mr^, T^) ® T3 ® (xj, X2)/A' 
+ E ('^')^3 ® n ® ix^, X.+i) ® (X- , X-+i)/A' 

i=l,3 

+ E i^)T3 ® T| ® (x., X.+i, X.+2) ® (X- , X-+1, X-+2)/A' 

i=5,8,ll 

+ zg;^($)T3 ® ® xi4 ® xu/A' + ^2) ® (Tf , T^) ® e ® evA^ 

+ ^($)T3 ® Tg^ ® e ® evAl 



(47) 



(4J 



In the same way, the Ti 2,3 scalar mass matrix can be written as 



[171. 



TJij 



0'ri2^3/2 



2 

3/2 

^T{1) + '^T22"^3/2 



3/2 



kpaiw? 



3/2 



A;t0.i^3/2 
A;tOi^3/2 

"^7^(2) + ^T33"^3/2/ 



(49) 



where fcy is a complex parameter whose magnitude is of order one, and it is the only new source of 
the CP violation in our model. The parameters Oipij are linear combinations of a^a^'s. In numerical 
analysis, we use the parameter Aai which is given by A 

In order to estimate the magnitude of FCNC phenomena, we move to the super-CKM basis by 
diagonalizing quark and lepton mass matrices including next-to-leading terms. For the left-handed 
down-type squark and slepton, we get 



{ml 



and for the right-handed down-type squark and slepton as 



l~2 \{SCKM) 



(50) 



(51) 



where the mixing matrices Vp and Up are given in Eq. ( ICllI) in Appendix C. 
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Let us study scalar trilinear couplings, i.e. the so called A-terms. The A-terms among left-handed 
an d rig ht-handed squarks (sleptons) and Higgs scalar fields are obtained in the gravity mediation 



as 



168 



hijLjRjHK = ^ hfpj^LjRjHK + h\^]^LjRjHK, 

K=5, 45 



(52) 



where 



''UK 



hfjiLjRiHK = -{yLjK)LjRiHKF''^K'^'^d^,Ku 
-{yiMK)LjRiHdF'''^K'^^d^,KMj 
-{yijK)LjRjHKF'^''K''^d^,KH^, 



(53) 



and Khi^ denotes the Kahler metric of Hx- In addition, effective Yukawa couplings of the down- type 
quark ijux are written as 



(54) 





/o 


09/^2 0^ 






/o 





o\ 


yijK = yi 





ag/ a/6 



















0) 


LR 












LR 



then we have 



h 



(y) 

UK 



m 

A 



F'^vVe 

\Q 0^ 







(0 





\ 


+ 













A 








LR 




^0 








(55) 



LR 



where F'^^ = F'^^/ai and F''- / A = 0{m^i2). 

By use of lowest level of the Kahler potential, we estimate h^^\ as 

hfjl = yuK{Af + A^j), 



(56) 



where we assume = A2 = A^ = F"'/(ajA) ^ 0(777,3/2). The magnitudes of Af = A2 and Af are 
also 0(7713/2) • Furthermore, we should take into account next-to-leading terms of the Kahler potential 
including Xi- These correction terms appear all entries so that their magnitudes are suppressed in 
0{a) compared with the leading term. Then, we obtain 



' ^LRll^'d Ci 2 ^LRlZ^d 



V^Lsi^rf «Li?32^rf C2mbir)J 



7.2 



(57) 



LR 



where al^jj are linear combinations of a^aes, and Ci and C2 are of order one parameters. Moving to 
the super-CKM basis, we have 



(-2 -siSCKM) _ ttU2 n t/ ^ _ , 



fO{d^v,) 0{~a'v,) 0{a^vA 



(58) 
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Similarly, for the charged lepton, 

/0(a%) 0{a''vd) 0{a\d)\ 



(59) 



C. Renormalization group effect 



In the framework of the supergravity, soft masses for all scalar particles have the common scale 
denoted by msusY, and gauginos also have the common scale mi/2. Therefore, at the GUT scale 
mcuT, we take 

Mi(mGUT) = M2{mGVT) = M3(mGUT) = mi/2 . (60) 
Effects of the renormalization group running lead at the scale mw to following masses for gauginos, 

M,{mw) ^ ^li!I^M,(mGUT). (61) 



Taking into account the renormalization group effect 169l | on the average mass scale in m 



m 



and with neglecting Yukawa couplings, we have 



ml^i^w) - mlimcvT) + 0.5M2(mGUT) + OMM^{mGVT) ^ "^iusY + OMml/^, 
^Inif^w) ^ "^r("^gut) + 0.15Mi^(mGUT) ^ "^iusY + 0.15mi/2 , 

ml^{mw) ^ m^(mGUT) + 0.004M2(mGUT) + 0.4M2(mGUT) + 3.6M32(mGUT) (62) 
-^susY + 4.1m?/2 , 

- "^/{("^gut) + 0.015Mi (ttigut) + 3.6M|(mGUT) - "^susy + S-Tm^/a ■ 
For Yukawa couplings, the h — r unification is realized at the leading order in our model, however, 
the h—T unification is deviated when we include the next-to- leading order mass operators due to terms 



including H/^^, see Ref. 147l | for the detail. In that paper, we have calculated the renormalization 
group equations and observed fermion masses at the weak scale can be obtained when tan /3 is larger 
than two. Hereafter, we take tan /3 = 3 on the numerical analysis. 



IV. NUMERICAL ANALYSIS 



In this section, we perform numerical analysis to show that the 5*4 flavor model presented in the 
previous section can explain the like-sign charge asymmetry in the Bg — Bg system. In order for this 
calculation, we first define the MI parameters for down-type squarks 5^^, 5^'^, S^^^, and and for 
sleptons S^^, Sf^, 5f^ and, as 

















-[ 



\{SC'KM) 



) (SCKM) 
\ (SCKM) 



YSCKM)\ 
ySCKM) I 

J (SCKM) \ 
^^2^jiSCKM)j 



(mH 



- diag(m^ 
diag(m? 



(63) 
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where rriq and are average squark and slepton masses with the values given below. 
In the numerical analysis, we fix the following parameters 



m3/2 = 430 GeV, = 880 GeV, = 520 GeV, 
Ml = 135 GeV, Ma = 270 GeV, M3 = = 1 TeV, (64) 

which are derived from the universal relation at the GUT scale 

"^i/2(mGUT) = '^3/2("^gut) = "^susY ("^gut) = 430 GeV, (65) 

by through the renormalization group effect discussed in the section III-C. This universal value is 
taken to be consistent with t 
Atlas Collaboration of LHC 



le lower bound of the gluino mass, which has been reported recently at 



159 



16 1| . For the other parameters, we assume the following regions: 



/X = [500,1000] GeV, AaL = [-0.5,5], |Mi| = [0,2], arg(A;Tai) = [-vr,7r], 
as = [0,0.001], dTi2 = [0,0.1], dLmj = [0,0.01], (66) 



with tan/3 = 3, ci,2 = 1, and kp = 1. In Eq. (]66l) . the number of left-handed and right-handed sides 
in braces denote the minimal and maximal values, respectively. In our calculation, we neglect the 
diagonal elements of scalar masses dpi^i^2,3) and aT(ii, 22,33)- The leading contribution to the parameters 
a in the soft-terms are ai as a^Kj ~ a/oiOs and dxu — cli- As given in Appendix [Dl the SUSY 



contribution to M^^^^^ is estimated as 



-9.4 + (5,^")y + 7.9(5,^^)23(5^)23}, (67) 

for X = rn^/m^- ^ 1.3, and similar for mixing. The coefficients in front of MI parameters for 
K meson mixing are —0.59,554, —183, 114, respectively. Since the LR terms are strictly 

constrained by 6 —t- 57 as seen in Appendix |Dl the terms gives larger contribution to 

Ml2^^^ . Among them, since ~ kpa\ < 10"^ in our parameter region given in Eq. fl66|) . the 

first term gives the dominant contributions. The approximation form of the LL parameters 

i^d^)ij given by 

2 



"^3/2 



V2 {ei^V^dk*Tai - ei^V^^kra^) , (68) 



i5^%s ^ -efs^aL + - ef,—^ kra, ^ ^y2KsfcTai, (69) 



.2 , /TT , /7r\ ^2 



{6^%s ^ -e',,AaL + ^ ( + e1,^-^] kra, ^ ^v^KdMi, (70) 



m? I 2 ^^2 / m'^q 
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where in the last approximation of each expression, we have neglected the first term proportional 
to ^^3 23 — 0.005 and Aa^, ~ 1. Notice that the MI parameters are expressed in terms of the CKM 
elements, and that both (5|'^)i3 and (5^^)23 have the same phase structure kxai, which is only the 
new source of the CP violation in our model. These are the typical feature of our 5*4 flavor model. 
By using these expressions, one finds that the — mixing induced by (5|'^)i2 is more suppressed 
by additional factor 9fj. 

The cEDM for the strange quark is estimated from the formula in Appendix |D] as 

erff ~ 10-'''lm[{S^%,{S^%,{6D32] ecm ~ 10-'Hm[{6!l%s] ecm, (71) 

for X ~ 1.3, {S^3)lr ~ 10~^ and {S^2)rr ~ 10~®. Therefore, (523)^^ is not constrained by cEDM. As 
for the b —> S'y process, one can see that (523)^^ should be strongly suppressed while {622)11, rr have 
an additional suppression factor mb/nig ~ 10"^. In our numerical calculation, we take diRij < 0.01 so 
that 6 — )■ 57 is well suppressed, and the allowed region of |((5|'^)23| is also small enough as mentioned 
below. 

First we discuss the allowed regions of the parameters {hg, hd), {hg, (Ts), {hd, (Td) defined in Eg. ( ITT]) . 
In our model, the parameters hd^sG^^"'^'° are estimated as 

j^d,SUSY 

hdc'-^ = -^fd;sM- ^ (27 - ^25)(5,^^)?3, (72) 

j^s,SUSY 

= -^f^ ^ (1.7 + z0.06)(5i^)^3, (73) 
M12 

where the MI parameters {Sd^)ij reflect the flavor symmetry, while the factors (27 — 225) and (1.7 + 
i0.06) do not. The ratio of (27-i25)/(1.7 + i0.06) is related to the CKM elements as (27-i25)/(1.7 + 
i0.06) ~ M^f^^/M^f^'^ ~ {Vtl/V;^)^. We obtain the ratio of hd and hs as 



ha 
h. 





27 - i25 






2 




2 


Vus 


2 


|1.7 + i0.06| 


\Vud 


2 ~ 


Vtd 


2 


Vud 


2 ~ 



(74) 



Therefore, the fact that the region hd — hg is favored reflects the flavor structure of the S4 flavor 
model. The CP violation phase 0^ is given by 

0, ~arg[-(l + M2-^)], (75) 

with neglecting the SM contribution. The CP phase sin^^ is bounded as | sin 0^1 ^ hg for hg < 1, 
and has the negatively-maximal valu e sin 0s — —hs at cr^ ~ 120°. This corresponds to the best-fit 



value {hs,as) = (0.5,120°) of Eq.(fT6|)|l62| 



Fig.l shows the plot in the 0^ — plane. The horizontal and vertical lines are the experimental 



values of one-dimensional likelihood analysis 



170 



0, = [-1.8,0.4] (rad), at 95% C.L., (76) 

and 2a range of in Eq.dH]), respectively. The blue (dark gray) and orange (light gray) regions 
denote 2a and 3a regions of A^i in Eq.dH]), respectively. By using Eq. (|T4]) and |1 + hsexp{2ias)\ ~ 
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FIG. 1: Allowed region in the c/i^ — A^^^ plane. The blue (dark gray) and orange (light gray) regions denote 



2(j and Sex regions of in Eq. ([8] 
experimental values of 2a region of 



The blue (black) error bars of the horizontal and vertical lines are 



and All. 
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Re[(5[;L)23l 



FIG. 2: Allowed region in the Re((^, 



LL\ 
d )23 



Im((5, 



LL\ 



J 23 



plane. The blue (dark gray) and orange (light gray) 



regions denote 2a and 2>a regions of A^, in Eq. ([8]). 



(1 + 0.5exp[2il20°]) ^ 0.8, we obtain —ali > 10 ^, and similar for af^. As a consequence we obtain 
the like-sign charge asymmetry as —A^i > 10~^, which is within 2a range of A^i of Eq.(|8]). 



Fig.2 shows the allowed region in the Re((5, 



)23 



Im((5^'^)23 plane. One finds from Eq.( 17^ that 



in order to obtain the best fit value {hg, ag) = (0.5, 120° 
must be opposite from each other, with Re((5^'^)23 — ±0.3 and lm.{6^^ 
region |(5^^ 



the sign of Re((5^^)23 and Im((5^ ;23 
=1=0.4. This allowed 



d J23 — =i=u.o aiiu imi^Orf ;23 

23 1 ^0.5 is small enough to suppress b ^ sj. The similar figure is drawn in the 
LL\_ -^rifin IT/ ./T/ I ~ Q22 times smaller area. 



Re(5^^)i3 - Im(5i^)i3 plane with \Vtd/Vts 



In Fig. 3, we predict t 
. The SM prediction 

;i - 5) X 10^= 



le difference a*; 



1511 « 



s,SM 
'■si 



— a 



^d,SM 
si 



a'^i, which will be measured at LHCb, as a function of 
(4.3 ±0.7) X lO"'* is also shown. We predict that 



— a'^i ~ (1 — 5) X 10^"^ in the 2a region of A^^. This will be a good test for our S4 flavor model. 
Since our model is based on the SU{5) GUT, above contributions in the quark sector affect 
;o the lep ton sector. Therefore, sleptons contribute to the LFV processe s and ED M of the electron 
17ll . ll72| . in which the experimental measurements give the upper bounds flTsl - llTsi . The Fig. 4 shows 
the relation of BR(/i — )■ 67) and the electron EDM. Within the MI parameters, {S^^)ij ~ {^d^)ji 
except for = (1,2), (2, 1) are relatively large in our model. Therefore one finds from Appendix 
E that the {6^^)21 term in which is enhanced by Mi/m^ ~ 10^, mainly contributes to ;U — )■ 67 
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process. As for the electron EDM, the terms with one small MI parameters dominates. The largest 
contributions are approximately estimated as 



BR(/i ^ 67) 



|4/e| 



3 X 10 

ai Ml 
An mj 



- 



An mj ytri^ 
11 /520GeVy 



2/2n(a;i 



Ml y/ |(^f^)2i| 
V mi J Vl35GeVy V IQ-^ 



Im[(5r)i3(5f^)3i]/3.(x 



+ Im[(5,^«)i2(5r)23(5r)31 + (5r)l3(5r)33(5f^)3l]/4„(Xl) 



^RR 



LR\ 



^RR\ 
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10-5 J \ 0.1 



(77) 



+ ■ ■ 



(78) 



The value of is of order m^i2Vd/rnj x a\j^^j < 10 ^ . Therefore we find that BR(/i — 67) and 

electron EDM can be close to the present experimental bound as shown in the figure. 
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The b — s transition by (5^^)23 in the quark sector simultaneously induce the LFV r decay r fij 
by (5^^)32. The dominant contribution is estimated from Appendix E as 



BR(r -> ^7) ~ BR(r ^ fii^^i^^)-^ l^— 



-/iMi tan /3 



2 



/520GeV\V Ml V Aan/^V / \2 



and similar for r — )■ 67 decay. Therefore for large term, r — )■ /i7 can be close to present upper bound 
given in Eq.f l2T]) . By using the expression Eqs. fl69|l and flTOj) . we obtain the relation of BR(r — > yU7) 
and BR(r — )■ 67) depending on the Cabibbo angle Ac — 0.22 as follows: 



ML^.M!ke.MM.^.A^. 05 (80) 
BR(r ^ /i7) - |(«32P " 1(5^^)23^ " " ^ " 

Therefore we conclude that there exist the parameter region which can explain the like-sign dimuon 
asymmetry A\i in the 6*4 flavor model, and in this case we predict that the LFV r — > yU7 decay can be 
so large that future experiments will reach, and the ratio of LFV of r decays, r ^ 67 and r — )• /i7, 
depends on the Cabibbo angle Ac. 



V. SUMMARY 

Recently the D0 Collaboration reported the like-sign dimuon charge asymmetry in hh — ?■ 
H^H^X decay processes. Their result shows 3.2a deviation from the standard model prediction. One 
promising interpretation of this result is that there exist additional contribution of new physics to 
the CP violation in Bg — Bg mixing process. In the effective Hamiltonian of the neutral Bg meson 
system, there are three physical quantities iF^gl) \^i2\ the CP phase 0^ = arg(— Mf2/ri2)- 
In order to obtain large CP asymmetry in the neutral Bg meson system, additional contributions 
from new physics to at least one of these three quantities are required. Within these possibilities, 
one can consider new physics that the absorptive part can be enhanced. However in general 
supersymmetric models, the gluino-squark box diagrams give the dominant contributions to Bg — Bg 
mixing, which do not affect |F^2l- Therefore in those models, new physics contributes to |Mf2| and 

In this paper we have considered an SU (5) SUSY GUT with 6*4 flavor symmetry. In this model, 
the Cabibbo angle, Ac ~ sin 15°, of the quark sector is given by a difference of 45° from up sector 
and 60° from down sector due to the Clebsch-Gordan coefficients at the leading order. As for the 
lepton sector, the tri-bimaximal form is generated in neutrino sector. These are consequences of the 
5*4 flavor symmetry. Since the matter multiplet T(10) and Fiji) are embedded into 2 + 1 and 3 of the 
5*4 group, respectively, the scalar masses of right-handed down- type squark and left-handed slepton 
are degenerated at the leading order, while those of Ti^2,3 fields are degenerated in the first two 
generations. Moreover for scalar mass matrix of Ti 2,3 fields, the relation (m|.)i3 = (m|,)23 oc kxai 
holds due to the 5*4 symmetry. The factor kxCii in the scalar mass matrix is assumed to be the only 
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additional complex parameter in our model, which is responsible for the CP violation in the neutral 
Bg meson system via gluino-squark box diagrams. As a consequence, the mass-insertion parameters 
(^i^)i3 and {6^^)23 have approximately the structure of VugkrCii and Vudkro-i, respectively. 

We have shown that the like-sign charge asymmetry A^i is in the 2a range of the combined result 
of D0 and CDF measurements. Since the relation between two CP phases sin (pa — sin (ps holds due 
to 5*4 flavor symmetry, and it can be large, we obtain large wrong-sign and like-sign asymmetry: 
{afi^l ~ {A^il ~ 10~^. The SUSY contributions in the quark sector affect to the lepton sector because 
of the SU{5) GUT relation ~ {^e^)ji- In the parameter region allowed by A^^, we have 

two predictions in the leptonic processes: (i) Both BR(/i — )■ 67) and the electron EDM are close to 
the present upper bound. Therefore, the MEG experiment 173| will be a good test of our model, 
(ii) The LEV r decays, r fij and r — )• 67, are related to each other via the Cabibbo angle Ac: 
BR(r — 7- e7)/BR(r — )■ fi'j) ^ A^. This is also testable at future experiments such as superKEKB. 
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Appendix A: Multiplication rule of S4 



The 5*4 group has 24 distinct elements and irreducible representations 1, 1', 2, 3, and 3'. All of 
the 5*4 elements are written by products of the generators bi and d^, which satisfy 



(61)^ = {d^Y = e, 6/4(61)^^4 = 61, d^bidi = 61(^4)^61 . 
These generators are represented on 2, 3 and 3' as follows. 



bi 



bi 



bi 



uj 


/o 1 \ 

10 , 

/o 1 \ 
1 



di 



1 

1 



dA 



dA 




on 2, 



on 3, 



on 3'. 



(Al) 



(A2) 



(A3) 



(A4) 



The multiplication rule depends on the basis. We present the multiplication rule, which is used in 
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this paper: 



ai 
a2 

ai 
ai 

a2 

V«3/ 



hi 
S 

b2 
Ws 

(h\ 

b2 

W 3' 

AA 

b2 



(ai6i + 0262)1 © {-ctih + ct2&i)i' 



aib2 + 0261 
ai6i — 0262 



02^1 

-|(V3ai62 + ct2&2) 
y KVSaifeg - 02^3) J 

^ aibi ^ 

^{V3a2b2 - aib2) 
\-^{VSa2b3 + 0163)/ 

{aibi + 02^2 + 03^3)1 3 



/ CLlbl 

i (1/80262 - 0.1^2^ 



\-|(y3a2&3 + 01^3)/ 3, 

/ 02^1 \ 

-|(A/3aife2 + a.2&2) 
y Kv^aifeg - 0263) / 3, 

^(a2&2 - 03^3) 
^^(-2ai6i + 0262 + a3&3 



'^0263 + 03^2^ 




'^03^2 


- 02^3^ 


0163 + 0361 




0163 


- asbi 


\a1b2 + a2&i y 


3 


\a2bl 


- 0162/ 







(bA 


02 


© 


b2 




3' 








(bl\ 


a2 


© 


b2 




3 





{aibi + a2&2 + dsb 



'3)1 



^(0262 - asbs) 
^^(-2ai6i + 0262 + a3&3 



'^0263 + «3&2^ 




'^0362 


- 02^3^ 


0163 + 0361 


© 


0163 


- 0361 


\a1b2 + a2&i / 


3 


\a2b1 


- ai&2/ 



(ai6i + 0262 + a3&3)] 



^(2ai6i - 0262 - asbs 
^(0262 - 03^3) 



'^0362 


- ^2^3^ 




'^a2&3 + 0362^ 




- asbi 
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0163 + a^bi 


\a2&i 


- ai&2y 
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\a1b2 + 0261 / 



More details are shown in the review 27 1. 
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Appendix B: Next-to-leading order 



Parameters appeared in the down-type quark mass matrix with next-to-leading order 
These are exphcitly written as 



ei2 = -T^yAfcOsau^'d - 



aia5Vd, 



ei3 



I - l)yA„, + |(v^ + l)yA„2 1 oiois - ^yAfe 05014 



Vd 



^ai3ai4 



^21 = VAc-^aianVd, 

£22 = —yAidbduVd - 



^(V3-l)yAe,+^(V3+l)yA 



^23 = <! -^(^ - 1)Z/A„, + -^(^ + l)yA„2 1 oioia - ^VA.a^au 



+ 



eai = -yA^a^a^Vd + t/Aj^agOia^'d, 
eaa = yA^a^agVd- 



Appendix C: Lepton sector 



The mass matrix of charged lepton becomes 



/O -3yiXagV45/V2 \ 

-3yiAa9t'45/-\/6 
\0 y2ai2,Vd) 



then, masses are given as 



= , mj = elyiAagl^v^ , = |?/2ra?3''^d 



In the same way, the right-handed Majorana mass matrix of neutrinos is given by 








o\ 

y^X^K-y^a^K 
My 
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and the Dirac mass matrix of neutrinos is 



D 










-as/ 




(o 





o\ 






+ y2^u 














) 




\a5 


as 


0-5/ 



(C4) 



By using the seesaw mechanism Mi, = M]~,Mj^^ Md, the left-handed Majorana neutrino mass matrix 
is written as 



where 



+ a — |6 a — |6 ^ 

a — |6 a + |6 + |c a + |6 — 
ya — |6 a + |6 — |c a + |6 + |cy 



M ' ?/fA2A + 2/2^a4A' |/fA2A-2/fa4A' 

It gives the tri-bimaximal mixing matrix f/tri-bi and mass eigenvalues as follows: 

7-r _ _J_ J_ L 

tvtri-bi - V3 V2 ' 

\ L ^ J- 



(C5) 



(C6) 



3a 



c . 



(C7) 



The next-to-leading terms of the superpotential are important to predict the deviation from the 
tri-bimaximal mixing of leptons. The relevant superpotential in the charged lepton sector is given at 
the next-to-leading order as 

Awi = i/A„(Ti, T2) ® (Fi, F2, F3) ® ixi, X2) ® (Xii, X12, X13) ® H-^/A^ 

+ l/A.(Ti, T2) ® (Fi, F2, F3) ® ixi, X2) ® (X5, X6, X7) ® ^45/A' 

+ l/A,(Ti, T2) ® (Fi, F2, F3) ® (xii, X12, X13) ® Xl4 ® i^45/A' 

+ 1/A,T3 ® (Fi, F2, F3) ® (X5, X6, X7) ® (X8, X9, Xio) (^H-^® /A' 

+ VAfTs ® (Fi, F2, F3) ® (X8, X9, Xio) ® (Xii, X12, X13) ® ^45 ® /A' . (C8) 

By using this superpotential, we obtain the charged lepton mass matrix as 



' ^11 — 2 



M, 



\ 



ei3 

621 ^ + 622 £23 



(C9) 



where and rrir are given in Eq. ( ]C2I) and ej/s are given as relevant linear combinations of auais. 
The explicit forms of e^j's are given by replacing y/^J^ with — I/a^ in Cij, which are presented in 
Appendix B. The charged lepton is diagonalized by the left-handed mixing matrix Ue and the right- 
handed one Ve as 



VIm.Ue = Mf 



(CIO) 
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where M?'*^^ is a diagonal matrix. These mixing matrices can be written by 



/ cos 60° sin 60° o\ 



Up 



( 



-sin 60° cos 60° 

V 1/ 

/ 1 f a 

1 a 



1 

A2 " 



17 



a 



+ ^ -a - ^ 1 



A A 

Taking the next-to-leading order, the electron has non-zero mass, namely 



(Cll) 



(C12) 



Appendix D: Formulae for quark sector 



Here we will give formulae for quark sector which are used in our analysis. The SUSY contribution 
by gluino-squark box diagrar a to the d ispersive part of the effective Hamiltonian for M — M mixing 
(M = Bd, Bs) is given by [l66|, [l76 | 



M.SUSY 
12 



al 2 



216m? 3 



A/ 



+ {5^%} {24x/6(x) + 66/6(x)} 



384 



M 



rrij + rrii 



M. 



M 



rrij + rrii 



36 ) Ux) 



•' I yrrij + rrii 

2 



M, 



M 



rrij + rrii 



- 84 } Ux) 



(Dl) 



where x = m'^/m'^ and the loop functions are defined as 



6(1 + 32;) logo; + x^ - dx"^ - 9x + 17 
6(x- 1)5 ' 

6x(l + x) logx — x^ — 9x^ + 9a; + 1 
3(x - 1)5 ■ 



(D2) 
(D3) 



For M = K,Bfi,Bs meson system, the generation indices of down-type quarks correspond to 
(ij) = (1,2), (1,3), (2,3), respectively. 

For b ^ s'-f decay, the Branching Ratio (BR) is given by 



BR(6 ^ .7) = 



m,G3ix)iS^%, + mgG,ix)iS^'')23\' + (L ^ R) 



(D4) 
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where tb is the hfetime of the B meson, and the loop functions are defined as 

1 + 4x — Sx^ + Ax loff X + log x 



Gi(x) 
Gsix) 



2(x-l)4 

-1 + 9x + — 17a;^ + 18a;^logx + Qx^ log a; 

12(a; - 1)5 ■ 



(D5) 
(D6) 



The chromo EDM of the strange quark is given by 163 1 



in rn?~ \ 3 
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(D7) 



where c is the QCD correction. We take c = 0.9. The functions A'^i(x) and N2{x) are given as follows: 

3 + 44a; - 36x2 - I2x^ + + 12a;(2 + 3x) log x 



Ni{x) 



6(x- 1)6 

10 + 9a; - 18a;2 - + 3(1 + 6x + 3^^) log x 
3(x- 1)6 ■ 



(D8) 
(D9) 



Appendix E: /x — )- 67, r — ;> 67 and t ^ iJ.j 

In the framework of SUSY, LFV effects originate from misalignment between fermion and sfermion 
mass eigenstates. Once non-vanishing off-diagonal elements of the slepton mass matrices are gen- 
erated in the super-CKM basis, LFV rare decays like ii — )■ ij'y are naturally induced by one-loop 
diagrams with the exchange of gauginos and slept ons. The decay ij — »• ij'j is described by the dipole 
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176 



178| 



operator and the corresponding amplitude reads [171, 

T = mi^e%{p - q)[iq''axu{ALPL + A^Pij)]ui(p) 



(El) 



where p and q are momenta of the initial lepton ii and of the photon, respectively, and Al^fi are the 
two possible amplitudes in this process. The branching ratio of ii — ij'-f can be written as follows: 



BR{ii ijUiUj) G] 



{\At? + \A 



r\ ) 



In the mass insertion approximation, it is found that 166 1 



At ^ 



012 



4tt m?- 
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llM2 



L(M|-^2 



f2n{x2,Xf,) + f2c{x2,X^) 



+ tan^ 9\Y jjiMi 



fsniXi) f2n{Xl,X^ 
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a 



1 i^n. 



47r m? 



M 



— 



ai 
Air 



^^/.M,tan/3^^='"(^^) '^^"^^^'^ 



2 

m- 



2 f2n{Xi 



(E2) 
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where 9w is the weak mixing angle, xi,2 = ^i,2/''^|) = /^^/'"^l fi(c,n){.x,y) 
fi(c,n){y)- The loop functions /j's are given explicitly as follows: 



fi{c,n) (•^) 



f2n{x) 
fsnix) 



— + 4x + 1 + 2x(x + 2) logx 

4(1 ' 
1 + 9x — — + 6x(x + 1) logo; 
3(1 

-x^ - 4x + 5 + 2(2x + 1) logx 
2(1 ■ 



(E3) 



Appendix F: Electron electric dipole moment 

The mass insertion parameters also contribute to t he electron ED M through one-loop exchange 



of binos/sleptons. The corresponding EDM is given as 



166 



179 



180 | 



dp 



where k^l = 2,3, {6^^)33 = —mT-{AT + /i tan /3)/m|, and the loop function /4„ is given as 

-3 - 44x + 36x2 ^ ^2x^ - x^ - 12x(3x + 2) logx 



(Fl) 



6(1 -x)6 



(F2) 



Since components (i, 3) and (3, i) of are much larger compared to others in our model, domi- 
nant terms are given as 



— 2 S ^( — ai) fsnixi) + 0{ — (1 + 



/itan/3^_ _2 



rrii 



)aia ) f4n{xi) >. 



(F3) 
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